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Abstract 

This brief report studies the behavior of entropy in two recent models of cosmic evolution by J. A. 
S. Lima, S. Basilakos, and F. E. M. Costa [Phys. Rev. D 86, 103534 (2012)], and J. A. S. Lima, S. 
Basilakos, and J. Sola [arXiv:1209.2802j . Both start with an initial de Sitter expansion, go through 
the conventional radiation and matter dominated eras to be followed by a final and everlasting 
de Sitter expansion. In spite of their outward similarities (from the observational viewpoint they 
are arbitrary close to the conventional Lambda cold dark matter model), they deeply differ in the 
physics behind them. Our study reveals that in both cases the Universe approaches thermodynamic 
equilibrium in the last de Sitter era in the sense that the entropy of the apparent horizon plus that 
of matter and radiation inside it increases and is concave. Accordingly, they are consistent with 
thermodynamics. Cosmological models that do not approach equilibrium at the last phase of their 
evolution appear in conflict with the second law of thermodynamics. 
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I. INTRODUCTION 



As daily experience teaches us, macroscopic systems tend spontaneously to thermodynamic 
equilibrium. This constitutes the empirical basis of the second law of thermodynamics. 
The latter succinctly formalizes this by establishing that the entropy, S, of isolated systems 
never decreases, S' > 0, and that it is concave, S" < 0, at least in the last leg of ap- 
proaching equilibrium (see, e.g., The prime means derivative with respect the relevant 
variable. In our view, there is no apparent reason why this should not be applied to cosmic 
expansion. Before going any further we wish to remark that sometimes the second law 
is found formulated by stating just the above condition on S' but not on S". While this 
incomplete version of the law works well for many practical purposes it is not sufficient in 
general. Otherwise, one would witness systems with an always increasing entropy but never 
achieving equilibrium, something at stark variance with experience. 

In this paper we shall explore whether two recently proposed cosmological models, that 
start from an initial de Sitter expansion and go through the conventional radiation and 
matter eras to finally enter a never-ending de Sitter phase, present the right thermodynamic 
evolution of above. This is to say, we will see whether S" > at all times and S" < at 
the transition from the matter era to the final de Sitter one. By S we mean the entropy 
of the apparent horizon, Sh = ks A/{4:ipi), plus the entropy of the radiation, 5*^, and/or 
pressureless matter, Sm, inside it. As usual, A and ipi denote the area of the horizon and 
Planck's length, respectively. 

We shall consider the evolution of the entropy, first in the model of Lima, Basilakos and Costa 
[2! and then in the model of Lima, Basilakos and Sola (models I and II, respectively). 
Both assume a spatially fiat Friedmann-Robertson- Walker metric, avoid -by construction- 
the horizon problem and the initial singularity of the big bang cosmology, evolve between 
an initial and a final de Sitter expansions (the latter being everlasting), and from the obser- 
vational viewpoint they are very close to the conventional Lambda cold dark matter model. 

In spite of these marked coincidences the physics behind the models is deeply different. 
While model I rests on the production of particles induced by the gravitational field (and 

2 



dispenses altogether with dark energy), model II assumes dark energy in the form of a 
cosmological constant that in reality varies with the Hubble factor in a manner prescribed 
by quantum field theory. 

We shall focus on the transitions from the initial de Sitter expansion to the radiation dom- 
inated era, and from the matter era to the final de Sitter expansion. As is well known, in 
the radiation and matter eras, as well as in the transition from one to another, both 5" and 
S" are positive-definite quantities IJ]. 

As usual, a zero subscript attached to any quantity indicates that it is to be evaluated at 
the present time. 



II. THERMODYNAMIC ANALYSIS OF MODEL I 

In model I the present state of cosmic acceleration is achieved not by dark energy or as a 
result of modified gravity, but simply by the gravitationally induced production of particles 
2|. In this scenario the initial phase is a de Sitter expansion which, due to the creation of 
massless particles becomes unstable whence the Universe enters the conventional radiation 
dominated era. At this point -as demanded by conformal invariance jol- the production of 
these particles ceases whereby the radiation becomes subdominant and the Universe enters a 
stage dominated by pressureless matter (baryons and cold dark matter). Lastly, the negative 
creation pressure associated to the production of matter particles accelerates the expansion 
and ushers the Universe in a never-ending de Sitter era. The model is consistent with the 



A. Prom de Sitter to radiation dominated expansion 



The production of massless particles in the first de Sitter era induces a negative creation 
pressure related to the phenomeno logical rate of particle production, F^, given by Pc = 
— (1 + w)pTr/ (SH) -see [Sj for more general treatments of the subject. Here, p is the energy 
density of the fluid (radiation in this case), w its equation of state parameter, and H = a/a 



H+la + ^)H'(l-^]=0, (1) 



the Hubble expansion rate. As a consequence, the evolution of the latter is governed by 
cf. equation (7) in j^. 

Because the rate must strongly decline when the Universe enters the radiation dominated 
era, it may be modeled as Tr/lSH) = H/Hj with H < Hj, being Hj the initial de Sitter 
expansion rate. In consequence, for w = 1/3 (thermal radiation) last equation reduces to 

H + 2H' (^1-^^=0, (2) 

whose solution in terms of the scale factor reads 



with D a positive-definite integration constant. 



The area of the apparent horizon A = 47rf^, where = is the radius trivially 

reduces, in the case under consideration (a spatially flat universe), to the Hubble length, 
H^^. Accordingly, the entropy of the apparent horizon, Sh = ksTc/iipi H)"^, as the Universe 
transits from de Sitter, H = Hj, to a radiation dominated expansion is simply 

It is readily seen that Sh is a growing, 5*^ > 0, and convex, S'l > 0, function of the scale 
factor (the prime stands for d/da). 



In its turn, the evolution of the entropy of t 
determined with the help of Gibbs equation jl| 



le radiation fluid inside the horizon can be 



dS^ = d{ Y^Aj + P7 ^ ( Y ) ' (5) 



where 

2' 



P7 = Pl 



1 + A 



(6) 



Pi is the critical energy density of the initial de Sitter phase, = Dal, Vi = Pi/'^- 
its turn, a* denotes the scale factor at the transition from de Sitter to the beginning of the 
standard radiation epoch. 
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Likewise, the dependence of the radiation temperature on the scale factor is given by 

-1/2 



-cfr. Eq. (13) in 



a. 



1 + 



(7) 



In consequence, 



A-K piD 
i^o^ = — — a > . 



3 nr— 

Accordingly, + S'^ > 0, i.e., the total entropy -which encompasses the horizon entropy 
plus the entropy of the fluid in contact with it- does not decrease. In other words, the 
generalized second law (GSL), first formulated for black holes and their environment js] and 
later on extended to the case of cosmic horizons js], is satisfied. 



Let us now discern the sign of 5"^' + S''. As we have already seen, S'l^ > 0. As for 5"'', we 



insert = T/ (1 -|- Da^ 



-1/2 



into Eq.(|8]) and obtain 



Thus, 



in piD 



3 HjTi 



1 + 2Da^ 



> 



(9) 



(10) 



^1 + Da^y/\ 

Therefore, S"^' + S" > 0; that is to say, in the transition from the initial de Sitter expansion 
to radiation domination, the total entropy is a convex function of the scale factor. If it were 
concave, the Universe could have attained a state of thermodynamic equilibrium and would 
have not left it unless forced by some "external agent". The initial de Sitter expansion 
{H = Hj, and no particles) was a state of equilibrium, but only a metastable one for the 
Universe was obliged to leave it by the production of particles which acted as an external 
agent. 



B. Prom matter domination to the final de Sitter expansion 



The Hubble function of spatially fiat Lambda cold dark matter models obeys 



H 



(11) 
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where Hoo = constant denotes its asymptotic value at the far future. 

In a dust filled universe {w = 0) with production rate Tdm < 3iJ of pressureless matter the 
Hubble factor obeys 

3 / 1 ^dm 



^ ^2 w)=^- ^^^^ 

Comparison with the previous equation leads to Tdm/i'^H) = {Hoo/HY , i.e., Tdm oc . 
Thus, 

H^ = Hlp^a-^ + (Ik], (13) 
where VL^ = {Hoo/HqY = 1 — Qm = constant > 0. 

Recalling that Sh = keA/iAeli), it follows = -2TTkB H' /{iliH^) with 

H-^J-H„ '1-"^')°-: . (14) 
2 °|(1 - SJA)a-3 + ^aI"^ 

Then 

S', = 6^ , (l-^A)a-^ (15) 



As for the entropy of dust matter, it suffices to realize that every single particle con- 
tributes to the entropy inside the horizon by a constant bit, say ks- Then, Sm = 
ks^r^n, where the number density of dust particles obeys the conservation equation 



n' 



{n/{aH))[Tdm - 3H] < with T,^ = SH^QJH > 0. 



Thus, 



S' - —kn — 



dm 



3H 



3H' 



(16) 



Since Tdm — 3i/ < and H' < the sign of S'^ is undecided at this stage. To ascertain it 
consider the square parenthesis in ( !T6ll and multiply it by aH/3. One obtains 

Tdm — 3H 



aH 



- 3H' 



(17) 



In consequence, S'^ > and the GSL, S"^ + S*^ > 0, is satisfied also in this case. 



6 



Let us now consider the sign of S'l + S'^ in the hmit a — )■ oo. From 5*^ = — {H' /H^) 

pi 

it follows, 

S'^ = -^-^i-AHH" -3H-]. (18) 



^pi 

In virtue of f lT^ we get 

H H" = -H' I <^ - - ^A)a-' + n^]a~' + (3/2)(l - I^a)^ a'^ 

2 ° 1 (1 - nA)a-3 + ' ^ ^ 



whence 



2 " ° 1 (1 - nA)a-^ + 



(20) 



Thereby, in view of ( !T8|) we get 5*^' < 0. 



As for the sign of S'^ it suffices to recall, on the one hand, Eq. ( 1T6|) and realize that 
^dm{ci — )■ oo) = 3HqQ\/ Hoo = SHoo and that H'{a — )■ oo) 0; then, S'^(a — > oo) = 0. 
And, on the other hand, that S!^{a < oo) > 0. Taken together they imply that S'^ tends to 
zero from below, i.e., that ^^(a — ?■ oo) < 0. 

Altogether, when a — )■ oo one has S'/l + S'^ < 0, as expected. Put another way, in the 
phenomenological model of Lima et al. the Universe behaves as an ordinary macroscopic 



system 



10| : i.e., it eventually tends to thermodynamic equilibrium, in this case characterized 



by a never-ending de Sitter expansion era with Hoo = Hq\/ VLj^ < Hq. 



III. THERMODYNAMIC ANALYSIS OF MODEL II 



a 



The model of Ref. [3| is based on the assumption that in quantum field theory in curved 
spacetime the cosmological constant is a parameter that runs with the Hubble rate in a 



specified manner [ll|, ll2|. As in the previous model, the vacuum decays into radiation and 
nonrelativistic particles while the Universe expands from de Sitter to de Sitter through the 
intermediate eras of radiation and matter domination. However, as said above, the physics 
of both models differ drastically from one another. While model I dispenses altogether with 
dark energy, model II assumes dark energy in the form of cosmological term. A, that evolves 
with expansion. 
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According to model II, the running of A is given by the sum of even powers of the Hubble 
expansion rate, 

A{H) = Co + 3uH^ + 3a— 2 , (21) 

where Cq, a and ly are constants of the model. The absolute value of the latter is constrained 
by observation as |z/| ~ 10^'^. At early times the last term dominates, and at late times 
{H <^ Hi) it becomes negligible whereby (pTl) reduces to 

A{H) = Ao + 3uiH^ - Hi) (22) 

with Aq = Co + SuHq. 

At the early universe, integration of the field equations gives for the _Hubble function, the 
energy density of radiation and of vacuum, the following expressions 



H{a) = \l^ ^[ {^ = {l-u){l + w)), (23) 



(l-z/)2 Da'^P A l-z/i/Da^^ + l 

where D{> 0) is an integration constant. 



As in model I, at the transition from de Sitter to the radiation era, the first and second 
derivatives of Sh and (with respect to the scale factor) are all positive. Thus, at this 
transition the GSL is fulfilled but neither the radiation era nor the subsequent matter era 
correspond to equilibrium states since at them 5"' > 0. 

By integration of the field equations at late times it is seen that the transition between the 
stages of matter domination to the second (and final) de Sitter expansion is characterized 
by 

H{a) = J(l-fiAo)a-3(i-) + Qao-z/ , (25) 

Pm = Pmoa"'^'"'^ and pM = paq + P"^^ [a-^^^-"^ - l] , (26) 

where fl\o = 87rGpAo(3ifQ)~^. From the pair of equations ( l26l) we learn that dust particles 
are created out of the vacuum at the rate Tdm = ^H. 
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As in the previous model, 5*^ > and 5*^' < 0. However at variance with it, the matter 
entropy, Sm = ks^rj^n oc H^^n, decreases with expansion and is convex. This is so 
because, in this case, the rate of particle production, Tdm, goes down and cannot compensate 
for the rate of dilution caused by cosmic expansion. Nevertheless, as it can be easily checked, 
5*^ and 5*^' dominate over S'^ and S'^, respectively, as a — )■ oo. Thus, as in model I, the 
total entropy results a growing and concave function of the scale factor, at least at the far 
future stage. Hence, the Universe gets asymptotically closer and closer to thermodynamic 
equilibrium. 



IV. DISCUSSION AND CONCLUDING REMARKS 



The second law of thermodynamics constrains the evolution of macroscopic systems; thus 
far, all attempts to disprove it by means of "counterexamples" have failed. While it seems 
reasonable to expect it to be obeyed also by the Universe as a whole, a proof of this on 
first principles is still lacking. However, persuasive arguments based on the Hubble history, 
suggesting that, indeed, the Universe behaves as any ordinary macroscopic thermodynamic 
system (i.e., that it tends to a maximum entropy state), were recently given in 10|. On 
the other hand, in view of the strong connection between gravitation and thermodynamics 



-see e.g. 
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14l |- it would be shockiu 



that the Universe behaved otherwise. In this spirit 



we have considered models 21] and [3|, each of them covering the whole cosmic evolution 
(i.e., the two well-known eras of radiation and matter dominance sandwiched between an 
initial and a final de Sitter expansions), and consistent with recent observational data. In 
both models, the entropy, as a function of the scale factor, never decreases and is concave 
at least at the last stage of evolution, signaling that the Universe is finally approaching 
thermodynamic equilibrium. 



In principle, the initial de Sitter eras should be stable {H and 5* are constants when t — )■ — oo) 
but owing to particle production, which can be viewed as "external" agent acting on the 
otherwise isolated system, an instability sets in. Once the Universe gets separated from 
thermodynamic equilibrium it reacts trying to restore it -as ordinary systems do-, only that 
at lower energy scale. This is finally achieved at the last de Sitter expansion. Given that two 
de Sitter expansions cannot directly follow one another an intermediate phase (comprised 
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by the radiation and matter eras) is necessary in between. 



As is well known, irreversible particle production, as is the case in models I and II, implies 
generation of entropy (see e.g., [sl, [l5|); something rather natural because the new born 
particles necessarily increase the volume of the phase space. Our analysis takes this into 
account in an implicit and straightforward manner via the F rates of particle production. 
These quantities modify the corresponding expressions for the Hubble factor and hence S" 
and S" . For instance, setting F,, to zero in Eq. ([T]) (which would kill model I) leads to 
D = and therefore to S'^ = S" = (Eqs. and ffTUl) . respectively). Likewise if, in the 
same model, one sets Tdm to zero, then 5*^ (Eq. f|T6|) ) decreases. Analogous statements can 
be made about model II if the parameter u (that enters the corresponding F rates) is forced 
to vanish (again, this would kill the model). 



en^opvof black hole honzons generalizes to S,^., 



A _ 1 ]r, A 



When quantum corrections to Bekenstein-Hawking entropy law are taken into account, the 

plus higher order 

terms [l6|, [l7| . Assuming this also applies to the cosmic apparent horizon, one may wonder 
up to what extent this may modify our findings. The answer is that the modifications are 
negligible whereby our results are robust against quantum corrections to the Bekenstein- 
Hawking entropy. We illustrate this point by noting that the expression for 5*^ of model I 
in the transition from the initial de Sitter regime to radiation domination presents now the 
overall multiplying factor ^ 1 — 8n{i+Da'^)^ \ • ^^^^ expression the second term is negligible 
on account of the quantity i^i in the numerator. It is noteworthy that the imposing of the 

£2 fj2 

condition S'/^ > leads to g^(/^£)a2)2 < 1- When this inequality is evaluated in the limit 
a — )■ 0, the upper bound on the square of the initial Hubble factor, H] < Sn/ipi follows. 
Thus, the nice and convincing result, that the initial Hubble factor cannot be arbitrarily 
large (its square, not much larger than Planck's curvature) arises straightforwardly from the 
quantum corrected Bekenstein-Hawking entropy law. 



Likewise, a study of S'f^ of model II in the transition from de initial Sitter expansion 
to radiation domination leads, in the same limit of very small a, to the upper bound 
H] < 27ra/[(l — 'y)ipi], i.e., to essentially identical result on the maximum permissible value 
of Hj. It is remarkable that in spite of being models I and II so internally different, they 
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share this bound. 



We conclude that models I and II show consistency with thermodynamics, and that their 
overall behavior (in particular, the reason why they evolve precisely to de Sitter in the long 
run) can be most easily understood from the thermodynamic perspective. Further, these 
results remain valid also if quantum corrections to Bekenstein-Hawking entropy law are 
incorporated. 



It would be interesting to explore the possible connection of the second law when applied to 
expanding universes with the "cosmic no-hair conjecture" [l8|. Loosely speaking, the latter 
asserts that "all expanding-universe models with positive cosmological constant asymptot- 
ically approach the de Sitter solution" [l9|. There is an ample body of literature on this 



-see e.g.. 



20| and 2l| and references therein. In the light of the above we may venture to 



speculate that the said conjecture and the tendency to thermodynamic equilibrium at late 
times are closely interrelated. Nevertheless, this is by no means the last word as the subject 
calls for further study. 



Before closing, note that the particle production does not vanish in the long run (rrfm(a — 
oo) = 3Hoo, and Tdm{(i — > oo) = j-^ifox/^Ao ~ ^ in models I and II, respectively). Then, the 
question arises as to whether it will be strong enough to bring instability on the second (an 
in principle, final) de Sitter expansion. Our tentative answer is in the negative; the reason 
being that in both cases the expansions tend to strictly de Sitter {H = constant > 0). 
However, a definitive response requires far more consideration and it lies beyond the scope 
of this work. At any rate, if instability sets in again, one should expect that the whole story 
repeats itself anew though at a much lower energy. 
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